Orbital Effects of In-Plane Magnetic Fields Probed by 
Mesoscopic Conductance Fluctuations 
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We use the high sensitivity to magnetic flux of mesoscopic conductance fluctuations in large 
quantum dots to investigate changes in the two-dimensional electron dispersion caused by an in- 
plane magnetic field. In particular, changes in effective mass and the breaking of momentum reversal 
symmetry in the electron dispersion are extracted quantitatively from correlations of conductance 
fluctuations. New theory is presented, and good agreement between theory and experiment is found. 
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A simplified view of transport in a planar two- 
dimensional conductor, as formed for instance by a semi- 
conductor heterostructure, suggests that when only the 
lowest quantized subband is occupied, an in-plane mag- 
netic field couples only to the electron spin, allowing 
the influence of an applied magnetic field to be sepa- 
rated into spin and orbital parts. However, the emerging 
picture of quantum transport in parallel fields 0, 0, 01 
has turned out to be surprisingly rich, indicating that 
even modest parallel (i.e., in-plane) fields can have sig- 
nificant orbital coupling, break time-reversal symme- 
try, and generate mesoscopic conductance fluctuations 
with field-dependent correlations — even without spin- 
orbit coupling or occupation of higher subbands. 

In this Letter, we use the high sensitivity of meso- 
scopic conductance fluctuations (CF's) to magnetic flux 
and time-reversal symmetry (TRS) to examine in detail 
the orbital effects of an in-plane magnetic field, B\\, in 
a quasiballistic quantum dot formed in a GaAs/AlCaAs 
2D electron gas (2DEG). Quantitative comparison of ex- 
periment and theory developed here allows the effects of 
,B|| on the electron dispersion in a planar 2DEG, includ- 
ing an anisotropic effective mass and a breaking of TRS 
(in spatially asymmetric confinement potentials), to be 
distinguished using various correlation functions of CF's. 
Effects of nonplanarity of the 2DEG are also included 
in the theory, and have distinguishable signatures in the 
CF correlations. The significance of the present work is 
to demonstrate experimentally that the effects of an in- 
plane field go far beyond Zeeman coupling, but cannot be 
characterized in terms of simple flux threading through 
the finite thickness of a 2D electron layer. Also, this 
study shows that phase coherent CF's can be used as a 
sensitive quantitative tool, much as one uses a supercon- 
ducting quantum interference device (SQUID). 

Effects of parallel fields on quantum transport have 
been investigated in 2D systems, including metal films 
silicon MOSFET's 5], and GaAs/AlGaAs 2DEG's 
013, as well as in ballistic focusing geometries Q. Those 



results were mostly interpreted in terms of flux threading 
due to finite thickness in the confined direction, surface 
roughness, or deformation of the Fermi circle due to the 
field. Subband depopulation, decoupling of bilayer sys- 
tems and diamagnetic shifts caused by Bn have also been 
observed using cyclotron resonance techniques [SjLmag- 
netoresistance measurements | lfJl . and tunneling as 
well as optical spectroscopy [T^]- Related investigations 
based on quantum-dot weak localization |3( were only 
sensitive to the breaking of TRS by a parallel field. Here, 
by using the full CF correlations, the effect of B\\ on the 
full electron dispersion is investigated, and the various 
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FIG. 1: Conductance fluctuations 8g(B^,B±,V) = 
g(B\\, B±,V) — {g(B\\, Bj_,V))v (color scale) through an 8/im 2 
dot with one fully transmitting channel in each lead measured at 
0.3K with a) shape-gate voltage V = —515 mV, b) By = and 
c) B± = 0, obtained from 2D slices of the three dimensional CF 
pattern, as indicated in d). 



2 



contributions are distinguished. 

Two quantum dots, with areas A = 8/im 2 and 3/im 2 , 
made on the same wafer, were measured and showed sim- 
ilar behavior. Data from the 8 /im 2 dot (see Fig. 3, in- 
set) will be presented in detail. The dots are formed by 
lateral Cr-Au depletion gates defined by electron-beam 
lithography on the surface of GaAs/AlGaAs heterostruc- 
tures grown in the [001] direction. The 2DEG inter- 
face is 900 A below the wafer surface, comprising a 100 A 
GaAs cap layer and a 800 A Alo.34Gao.66As layer with a 
400 A Si doped layer set back 400 A from the 2DEG. An 
electron density of n = 2 x 10 15 m~ 2 and bulk mobility 
/i ~ 14 m 2 /Vs (cooled in the dark) gives a transport mean 
free path £ e ~ 1 /im. Note that the 8 /im 2 dot contains 
of order 10 4 electrons. 

Measurements were made in a 3 He cryostat at 0.3 K us- 
ing current bias of 1 nA at 338 Hz. Shape-distorting gates 
were used to obtain ensembles of statistically indepen- 
dent conductance measurements 13] while the point con- 
tacts were actively held at one fully transmitting mode 
each (N = 2). In order to apply tesla-scale Bu while 
maintaining sub-gauss control of B± , we mount the sam- 
ple with the 2DEG aligned to the axis of the primary 
solenoid (accurate to ~ 1°) and use an independent split- 
coil magnet attached to the cryostat to provide B± [l4|. 
The Hall effect measured in a separate GaAs Hall bar 
mounted next to the quantum dot, as well as the loca- 
tion of weak localization minima in transport through 
the dot itself (visible Bn < 2T) were used to determine 
the offset in B± (i.e. the residual sample tilt), which was 
then corrected by computer control of the two indepen- 
dent magnets. 

The raw data consist of measured dot conductance 
g(Bu, B±, V) as a function of shape-distorting gate volt- 
age V (inner loop of multiparameter sweeps), B±, and 
B|| (outer loop, swept from — 2.5T to +4T over ~ 20 h), 
giving 20 independent shape, 15 independent B± and 
about 10 independent Bu samples. Conductance fluctua- 
tions are found by subtracting the gate- voltage averaged 
conductance over the measured range: Sg(Bu, B±, V) = 
g(BpBx,V)-{g(B h B ± ,V)) v . 

Figure 1 shows 2D slices of conductance fluctuations 
in the full 3D space of Bn, B±, and V. Note that be- 
cause gate-voltage-averaged conductance is subtracted 
from the fluctuations, weak localization effects on (g) are 
not evident in Fig. 1. A principal result is already evident 
in Fig. 1: The horizontally elongated features around 
~ in Fig. 1(c) show qualitatively that CF's are 
less sensitive to Bn in the vicinity of \B«\ ~ 0, giving a 
larger correlation field near |£fii| ~ 0, than at larger par- 
allel fields. This elongation, demonstrating reduced flux 
sensitivity near \Bu | ~ 0, is consistent with the Bu de- 
pendent effective mass and momentum reversal symme- 
try breaking terms of our theory. Effects of nonplanarity 
alone would result in a Bu independent correlation field. 
A quantitative analysis is presented in Fig. 2. 

The 2D slices in Fig. 1 also illustrate the fundamen- 
tal symmetries of conductance with respect to magnetic 



fields Bu and B±: when Bu — 0, conductance is sym- 
metric under inversion of B±, g(B±) = g(—Bj_) (see 
Fig. 1(b)); when B± = 0, conductance is symmetric 
under inversion of Bn, g(Bu) = g{— Bn) (see Fig. 1(c)). 
When both Bu and B± are nonzero, the symmetry of con- 
ductance requires the reversal of both fields, g(Bn , B±) = 
g(—Bu,—B±) (see Fig. 1(a)) ^{|. The fact that a nonzero 
Bu breaks the symmetry g(B±) = g(—B±) is a simple 
qualitative demonstration that By breaks TRS A 
quantitative analysis of this effect is presented in Fig. 3. 

To quantify the correlations of the various parameters 
used to generate CF's — including in particular Bu - 
we define the normalized correlation functions, 



C V {B\\) 
C h ( B \\) 



(Sg(B h V)Sg(B h V + v)) 

(6g{B h B ± )8g{B h B ± + b ± )) 

(W||)> 
(SgjB^SgjB^ + b^ 

vW(WW)' 



(1) 

(2) 
(3) 



where (...) is shorthand for (. . .)v,b±i i- e -> averaging over 
both gate voltage and B±, with B± sufficiently large to 
fully break TRS throughout the measured range. 

Theoretical expressions for the correlation functions in 
Eqs. 1-3 can be found using the effective 2D Hamiltonian 



H 2C - 7^-P±rt B \\) +PlP( B \\) + «W> (4) 

for electrons confined to a plane perpendicular to z 0|. 
Here, p = —iTiV — - A 2 d, with rotA 2 £> = B±, is the 
2D momentum operator in the plane, with component 
p± = p- [Bu x l z ]/Bu perpendicular to Bn, and it(r) is the 
impurity and dot confining potential. The middle terms 
in H2D arise from pj_-dependent subband mixing: the 
7(-B||) term lifts rotational symmetry with an anisotropic 
mass enhancement but does not break TRS; the /3(B^ 
arises from the broken spatial inversion symmetry of the 
heterointerface and breaks TRS in presence of Bn. 

Extending the analysis of Ref. [l| to higher temper- 
ature yields a correlation between conductance fluctua- 



tions at Bn and at Bu 



C h {B\0 = 



6|| of the form 
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in the unitary ensemble, where t~ s ], = NA/h is the es- 
cape rate from the dot, with A = 2irfi 2 /m*A the mean 
level spacing of the corresponding closed dot (effective 
electron mass m* = 0.067m e ), tJ 1 is an additional es- 
cape rate due to orbital effects of Bu , as discussed below. 
The exponent a equals —1 in the high temperature limit 
kT 3> (Tit~ s I, fiT^ 1 , ez), applicable in the present ex- 
periment, and —2 in the low temperature limit, where 
ez = g^BB is the Zeeman energy, with g = —0.44 for 
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GaAs. The difference between the high and low tem- 
perature regimes is caused by the necessity to average 
the interference contributions coming from electrons at 
different energies. For parallel fields with ez < 3fcT, ap- 
propriate for the present measurements, the deviation of 
Eq. (5) from the full expression is negligible jlfij . 

The additional escape rate tJ 1 due to Bn is given by 



a ) by (B|| = 2.2 T) ^ bj (B„ = T) bl (B„ = T) V^B, = T) 
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where r = [im* /e is the elastic scattering time in a dif- 
fusive dot or the crossing time r — m*L/pp in a bal- 
listic device, where L is the diameter of the device and 
Pf = h(2nn) 1 / 2 is the Fermi momentum. The £ term de- 
scribes effects of nonplanarity, including interface rough- 
ness and dopant inhomogeneities, and also breaks TRS. 

Writing the functions 7(0;) and (3{x) in Eq. (6) in terms 
of scale factors 7 and (3 and normalized functions g{x) 
and f{x) as 7(2;) = ^g{x) and (3(x) = /3/(x), we find 
g(x) and f(x) from self-consistent simulations of the het- 
erostructure 01 an d treat 7 and (3 as fit parameters. 
Below ~ 2 T , the normalized functions are well approxi- 
mated by g(x) k, x 2 and f(x) « a; 3 (see Fig. 2(b), insets); 
however, the full functions are used for all comparison of 
theory and experiment. We note that 7 and (3 can also 
be obtained from the heterostructure simulations, giving 
values in reasonable agreement with those obtained from 
the fits. 

Figure 2(a) shows experimental correlation functions, 
Chu , Cb ± , and C v , for representative parallel fields, as in- 
dicated. The corresponding characteristic voltage v c and 
fields 6jj and 6j_ are shown in Fig. 2(b) as a function of 
B\\ , as obtained from the half width at half maximum 
(HWHM) values of the correlation functions, indicated 
by dashed lines in Fig. 2(a). It is evident from Fig. 2(b) 
that both b c ± and v c are independent of B\\ within the 
error bars, in agreement with theory and previous exper- 
iments (An alternative procedure, not shown, for 
extracting these same quantities from the slopes of log- 
power spectra of CF's yields similar values for v° and 6| 
that are again independent of B\\ , within error bars.) 

In contrast, the parallel field correlation length, 
shown in Fig. 2(b) decreases substantially from its zero- 
field value on a field scale of ~ IT. Good agreement with 
theory is found: the solid curve in Fig. 2(b) is the best- 
fit (described below) theoretical HWHM correlation field 
for Cb y (i?||) obtained from Eq. (5). This decrease is due 
to the 7 and (3 terms in Eq. (6) and cannot be accounted 
for with the £ term alone. 

Symmetries of conductance in parallel and perpendic- 
ular fields are investigated in Fig. 3. We define the cross- 
correlation functions 



C±±(B l{ ) = 



(6g(B h B ± )Sg(±B h ±B 1 _)) 
y /(8g*{B h B ± )){6g*{±B\ U ±B ± )y 




0.0 0.8 1.6 0.0 0.4 0.8 20 40 
bn (T) b ± (mT) 



v (mV) 



b) 



0.4 



0.0 





9 / 




60 




5 



IT 2 ] / ■ 


W by theory 


50 


0.6 


20 
"10 


, f 


f \ hi expt. 


40 <„ 









A V 




0.4 „ 




12 3 




30 % 


3 

H 








20 


0.2 




I 




10 













0.0 



(7) 



1 



FIG. 2: a) Correlation functions C(,,. at ZJ| = 0, 2.2 T (open cir- 
cles, squares) as well as C), ± (open diamonds) and C v (open tri- 
angles) at By = 0. Half width at half maximum values give the 
characteristic voltage v c and fields and 6|, shown in b) (solid 
circles) as a function of By . Markers in b) refer to corresponding 
curves in a). The solid curve shows the three-parameter theory. A 
typical error bar is indicated. Insets: g(x) and f(x) used for fits 
(see text) as obtained from numerical simulations (solid curves) as 
well as quadratic and cubic low field approx. (dashed curves). 



With this definition, the first (second) subscript index of 
C refers to whether B\\[B±) is reversed when comput- 
ing the correlation function. Correlations for total field 

reversal, C (i-e., both B\\ and B± inverted) remain 

near unity for all parallel fields, as expected from the full 
Landauer-Buttiker (Onsager) symmetry (see Fig. 3(a)). 

Deviations from a perfect correlation C = 1 are 

small, indicating that the confining potential of the dot 
did not drift significantly over periods of a day. Fig- 
ure 3(a) also shows conductance fluctuations at (B\\,B±) 

and (B\\,—B±) (C_| ) become uncorrelated (C_| ~ 0) 

at parallel fields of a few tesla, indicating the field scale 
at which Bn breaks TRS. Within error bars, C |_ is in- 
distinguishable from C_| , as expected from Landauer- 
Buttiker symmetry, 5g(—Bu,B±) = Sg(Bn,—B±). The 
theoretical cross-correlation using Eq. (5), shown as a 
solid curve in Fig. 3(a), is in very good agreement with 
experiment data. 

Finally, we discuss the full correlation, Cj^{B\\), of 
CF's at B\\ with CF's at Bn + bn. Representative curves 
for Bn = 0, 0.6, 1, 1.4 T are shown in the lower part of 
Fig. 3 as a function of Bn + bn , along with best-fit the- 
ory curves based on Eq. (5). Besides the perfect cor- 
relation at 6|| = (Cfc ||= o = 1), there is an "echo" of 
correlations, both in experiment and theory, that occurs 
at fop ~ — 2-Bii. Within the present theory, this field- 
reversed correlation "echo" , C-2B« , is suppressed from 



4 



CD 



o 



1.0 
0.5 
0.0 


. " C— (B||) 

I ssiS*** 


f - C- + (B n ) 
vC + -(B„) 


— theory 
--■ LB invar; 




-2 -1 


1 
B„ (T) 


2 


1.0 
0.5 


- B„=0T 

I 




• exp. 
- theory 

- 


0.0 








1.0 
0.5 


B|| = 0.6 T 

I 




• GXp. 

— theory 


0.0 








1.0 
0.5 


- B„ = 1 T 


' ' ' — v-s 


• exp. 
& - theory 










0.0 








1.0 








0.5 


B|| =1.4T j 


Wr y 


\- theory 


0.0 










-2 -1 


1 
B || + b„ (T) 


2 



FIG. 3: Top: Cross-correlations of CF's at (B\\,B±) with CF's 

at (— Bn,— Bj_) (C , crosses), (— Bn,B±) (C_-|_, squares) and 

(Bu,—B±) (C-f — , triangles) as a function of Bn . Landauer-Biittiker 

(LB) symmetry for full field reversal gives C = 1 (Dashed line). 

Below: Correlation functions Cf, (By) at By = 0, 0.6, 1, 1.4T. 
Solid curves show theory based on Eq. (5). Typical error bars as 
indicated. 



unity only to the extent that parallel field breaks TRS. 
The agreement between theory and experiment, includ- 
ing the unusual asymmetric curves in Fig. 3, is quite 
good. A single, consistent set of three parameters (7, $, 



Q have been obtained from fits of Eq. (5) to 131 curves 
like those in Fig. 3, ranging over — 2.5T < Bn < 4T. 
We emphasize that all theory curves shown in Fig. 2 and 
Fig. 3 used this single set of three fit parameters and 
were not individually fit. The values obtained in this 
way were 7 = 11 ± 2 x 1CT 4 [m*]- 1 ^ 2 , j3 = 4 ± 4 x 
10- 4 [m*pF] _1 T- 3 , and C = 44 ± 8 x Hr 3 ^] -1 ? 1 - 1 , 
consistent within the error bars with values extracted for 
the 3 /1m 2 dot. We note that parameters {3 and £ obtained 
from the parallel-field-induced crossover from the orthog- 
onal to the unitary ensemble |3j are in good agreement. 
The self-consistent simulations give theoretical values of 
7 = 35 x 1(T 4 [to*]- 1 ! 1 - 2 , = 3 x 10- 4 [m*p F ] _1 T- 3 . 

In summary, orbital effects of an in-plane magnetic 
field Bn were experimentally investigated using the high 
sensitivity of CF's to magnetic flux in a large quasibal- 
listic quantum dot. Detailed quantitative comparison of 
correlations of CF's induced by Bn with theory devel- 
oped here reveal the mechanisms of coupling, including 
an induced anisotropic effective mass, the breaking of 
time-reversal symmetry due to the heterostructure asym- 
metry and effects of nonplanarity. In the present exper- 
iment, spin-orbit coupling is weak. On the other hand, 
the combined influence of stronger spin-orbit coupling Q 
and parallel fields is expected to yield interesting addi- 
tional features in the correlations and symmetries of CF's 
[l8|]. These remain to be investigated experimentally. 
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